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THE  FLOW  FIELD  BEHIND  A  SPHERICAL  DETONATION  IN  TNT  USING  THE  LANDAU - 
STANYUKOVICH  EQUATION  OF  STATE  FOR  DETONATION  PRODUCTS 


Prepared  by: 
M.  Lutzky 


ABSTRACT:  Calculations  have  been  made  of  the  flow  field  in  the  isentropic 
region  behind  a  detonation  wave  in  TNT,  using  the  Landau-Stanyukovich 
equation  of  state  for  the  detonation  products  (as  described  by  Zeldovich 
and  Kompaneets).  Adjustable  constants  in  this  equation  have  been  evalu¬ 
ated  by  imposing  ideal  gas  behavior  on  the  detonation  products  in  the 
large  expansion  (low  density)  limit,  and  by  fitting  to  an  experimental 
curve  of  detonation  velocity  versus  loading  density.  Calculated  values 
of  Chapman -Jouguet  variables  correspond  fairly  well  with  experimental 
values  at  various  loading  densities,  with  the  exception  of  the  tempera¬ 
tures,  which  seem  to  be  far  too  low.  This  is  connected  with  the  fact  that 
the  theory  predicts  an  upper  limit  to  the  loading  density  at  which  an 
explosive  will  detonate;  at  this  point  the  thermal  energy  vanishes  and 
only  the  elastic  energy  contributes  to  the  energy  of  detonation. 
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Calculations  of  the  airahock  motion  produced  by  a  spherical  TNT  explo¬ 
sion,  with  the  reaction  products  considered  to  be  gaseous,  have  given 
satisfactory  agreement  with  experimental  results  However,  the 
experimental  motion  of  the  explosive  interface  and  the  second  shock 
have  not  agreed  with  the  theoretical  calculations .  An  attempt  to 
clarify  these  discrepancies  has  led  to  consideration  of  the  Landau- 
Stanyukovich  solid  state  model  for  the  reaction  products  of  a  condensed 
explosive.  The  Iandau-ot&nyukovich  equation  of  state  has  been  utilized 
to  calculate  the  flow  field  in  the  reaction  products  behind  the 
Chapman -Jouguet  zone  -  the  so-called  Taylor  Wave  distribution  -  and 
the  results  are  presented  in  this  report  Preliminary  determinations 
of  this  distribution  have  already  been  used  as  initial  conditions  for 
the  calculation  of  the  subsequent  explosion  motion,  and  have  been 
reported  elsewhere. 
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INTRODUCTION 

The  theory  of  the  detonation  process  for  a  Y-lav  gas,  whose 
deto’atior  product  is  also  a  y-law  gas,  has  been  quite  completely 
worked  out  and  is  available  i>  many  places.1’3  In  particular,  the 
conditions  at  the  Chapman  -  Jougue+  state  can  be  derived,  a^d  it  can 
be  shown  that  the  detonation  velocity  i3  a  function  only  of  the  heat 
of  detonation  and  the  Y  for  the  detonation  products,  for  sufficiently 
large  detonation  pressures.  In  addition,  differential  equations  have 
been  derived  for  the  flow  behind  the  detonation  wave,  and  have  been, 
solved  for  certain  explosives  and  geometries 3,4 

The  theory  of  the  detonation  of  a  solid  explosive,  on  the  other 
hand,  is  in  a  much  less  satisfactory  state.  Experiments5’8  have  shown 
that  the  detonation  velocity  of  a  solid  explosive  depends  on  the  initial 
density,  unlike  the  detonation  velocity  of  gaseous  detonations.  Further¬ 
more,  the  explosion  products  of  condensed  explosives  are  obtained  at 
pressures  of  the  order  of  mega bars ,  and  at  densities  approaching  2  grams/ 
cm1,  under  which  conditions  their  behavior  becomes  extremely  complex. 
Consequently,  various  attempts  have  been  made  to  find  an  equation  of 
state  for  the  explosion  products  by  treating  the  highly  compressed 
gas  as  a  solid. 

The  first  such  attempt  is  due  to  H.  Jones,7  who  developed  an 
equation  of  the  Gruneisen  type,  based  on  the  Einstein  model  of  a 
solid,  of  the  form  p  *  Ae“av-B+fRT,  where  a,  A,  B  and  f  are  constants. 

The  equation  of  state  which  we  consider  in  this  paper,  however,  was 
derived  by  Landau  and  Stanyukovich  ,*  »*  who  also  approached  the  problem 
by  drawing  an  analogy  between  the  suate  of  the  detonation  products 
of  a  condensed  explosive  and  the  crystal  lattice  of  the  solid  state. 

It  is  well  known  that  the  energy  of  a  solid  body  has  a  two-fold  origin: 
it  is  made  up  of  an  elastic  energy  arising  from  the  binding  forces 
between  the  atoms  and  molecules  and  a  thermal  energy  connected  with 
oscillations  of  the  atoms  or  molecules  about  their  positions  of  stable 
equilibrium.,  Landau  and  Stanyukovich  have  attempted  to  describe  the 
behavior  of  the  detonation  product  by  considering  it  as  a  solid  with 
the  property  that  the  elastic  energy  and  the  elastic  part  of  the  pres¬ 
sure  are  predominant.  The  theory  has  been  described  and  expanded  by 
Zeldovitch  and  Kompaneet9®  so  that  we  refer  to  it  as  the  I3ZK  theory. 

The  purpose  of  this  paper  is  to  make  some  computations  using  the  ISZK 
equations  of  state,  and,  in  particular,  to  calculate  the  flow  field 
behind  the  detonation  shock  in  a  condensed  explosive. 
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THE  LSZK  EQUATION  OF  STATE 


For  the  sake  of  completeness ,  we  present  here  a  description  of 
the  I£ZK  equation  of  state,  and  a  derivation  of  some  of  its  properties. 


The  LSZK  equation  of  state  may  be  written8 
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where  P  ■  pressure 

E  *  energy  density  (per  unit  mass) 
v  *  specific  volume 

T  *  temperature 


sind  B,  Cyj  ,  Cv  aod  y  are  constants,  y  is  a  dimensionless  constant 
serving  as  a  polytropic  index  connected  with  the  intermolecular 
forces;  Cv  is  the  specific  heat  at  constant  volume;  Cyj  is  a  specific 
heat  associated  with  the  appropriate  lattice  vibrations;  and  B  is 
a  constant  having  the  units  ,  gm  calories  .  The  elastic 
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Eliminating  T  between  (l)  and  (2) ,  we  obtain  the  expression 


P  '  tr^  *  7T  l1  ‘  STV^T  }  • 

(3) 
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00 

O'  is  a  convenient  variable  which  will  be  used  in  this  report, 
terms  of  or,  (1)  and  (2)  may  be  written: 

In 
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(6) 


Another  convenient  parameter  which  we  will  find  useful  is  the 
quantity  y,  defined  as  the  ratio  of  the  thermal  part  of  the  pressure 
to  the  elastic  part: 
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(CyT/»v) 

(B/vV) 
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Clearly,  (5)  and  (6) 


may  now  be  written 


in  vhe  form: 


p  -  ~  (i  +  y) 
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ISENTROPIC  PROCESSES 

It  is  possible  tr  obtain  an  expression  for  the  pressure  of  the 
form  P  -  P  (p),  valid  for  iser.tropic  processes  of  an  LSZK  substance, 
by  combining  equation  (3)  with 

dE  -  dp  ,  (10) 

which  is  the  differential  equation  of  ar  isentropic  process  Dif¬ 
ferentiating  (3) ,  we  obtain 

dP  »  ~  (pdE  +  Edo)  +  Fv,pv"1  |i  -  ^ny  }  ^  (n) 

Using  (10)  to  eliminate  dE  from  (ll),  we  obtain 

dP  -  i  (|dp  +  Edp)  +  BYp^”1  (l  -  }  dp  (12) 

Solving  (3)  for  E  and  substituting  into  (12),  we  obtain  the  differential 
equation 


dP 

dp 


'1  +  0\ 
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which  has  the  solution: 


1+3 

P(p)  -  Kp"7-  +  BpY 


(13) 


(14) 


where  K  is  a  constant  of  integration.:  We  ar**  now  in  a  position  to 
obtain  expressions  for  E,  the  sound  speed  c  ,  the  temperature  T, 
etc-  as  functions  of  density  alone,  valid  for  isentropic  processes 
of  an  LSZK  substance.;  Thus,  putting  (14)  into  (3),  and  solving 
for  E,  we  obtain: 
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Similarly , 
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CHAPMAN - JOUGUET  CONDITIONS 

We  now  obtain  the  initial  conditions  at  an  LSZK  detonation,  in 
terns  of  the  variable  y.  We  consider  that  the  detonation  wave  consists 
of  a  shock  traveling  at  speed  D,  followed  iranediately  by  a  region  of 
isentropic  expansion.  The  region  of  chemical  reaction  behind  the 
shock  is  considered  to  be  infinitely  thin.:  Values  of  the  hydro¬ 
dynamic  parameters  in  the  undetonated  explosive  ahead  of  the  shock 
are  given  a  subscript  o  so  that  v0  is  the  specific  volume  of  the 

original  explosive.  We  first  obtain  ^  as  a  function  of  y.  From 

vo 

the  Rankine-Hugoniot  (R-H)  relations  at  the  shock,  we  have 


3l.JL  • 

v  D-u  ' 

where  u  *  particle  velocity 

D  ■  detonation  velocity. 

Using  the  detonation  property  D  ■  u  ♦  c 
•fuation  (l8)  becomes 

-  Ji  +  i 

v  c 

Another  R-H  relation  yields 

P  -  Du/v0 

which  may  be  written  in  the  form 


(18) 


(19) 


(20) 


(21) 


(22) 


Using  P-Bd^  (1  +  y)  and  the  equation  (20),  (22)  may  be  put  into  the 
form 


(23) 
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ca  may  be  obtained  as  a  function  of  p  and  y  by  eliminating  K  between 
(14)  and  (IT),  and  then  using  P  *  BpV  (l  ty)  to  eliminate  P;  the 
result  is 


c3  -  Bpv_1  (y  ♦ 


(24) 


Inserting  (24)  into  (23)  we  obtain  as  a  function  of  y: 


ifl. 


i  ♦  U  t-d _ 


(25) 


The  next  step  is  to  obtain  v0  itself  as  a  function  of  y;  thi6  will 
enable  us  to  solve  for  the  parameter  y  as  a  function  of  the  known 
quantity  v0  . 

We  write  the  R-H  equation  for  the  energy  in  the  form 


4  +  5  p  (v0  -  v)  -  a  +  5  Pv  (£■  -  1)  , 


(26) 


where  Q  is  the  chemical  energy  released  by  each  gram  of  explosive; 
and  using  (25),  this  becomes: 


E  -  Q  + 


§PV(  )  . 

'  Y  + 


(27) 


Eliminating  P  by  using  P  ■  Bp^  (l  +  y)  we  obtain 


E 
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i  +  2  ji_Lyll  (, 

Q  r\ _ V”1  V 


2vy_1  V  + 


TlUJ 

l  cr  J 


(26) 


Another  expression  for  e/Q  may  be  obtained  from  equation  (9).: 


e  .  b  + 

3  Q  VY-1  ^ 


(y-l)^ 


(29) 
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Equating  (28)  and  (29),  and  solving  for  v,  we  ODtai: 


.Y”^-  »?/_„,+  ^ 


v  ■  -  °  +  _±_  — C1  +  y).a_ —  -  . 

Q  *  V_1  v2(V  +  y  )y 


(30) 


Eliminating  v  between  (30)  and  (25),  we  obtain  t  e  expression: 


(5C-i  ;OT  +  J_  .  i2JLllirY-i  (x  *  Llj)  (n) 

V  l07  Y-l  2vT  /  li  W 


where 


w  =  v  + 


1  ♦  or 

■ — ;y 


(32) 


Since  v0  ,  the  specific  volume  of  the  solid  explosive,  is  a  known 
quantity,  we  may  solve  (31)  for  y,  by  an  iterative  process.  Since 
v  is  a  known  f unction  of  y,  by  virtue  of  (30)  ,  we  can  find  P  by 
using  the  expression 


(1  +  y); 


E  may  be  found  from  (28)  or 


T 

(29);  and  c3  may  be  found  from  equation  (24). 


The  particle  velocity  at  the  front,  u,  may  be  found  from  (20), 
and  is  given  by 


u 


c 


y 


(1  *  y) 


(33) 


Finally,  the  detonation  velocity  can  be  found  from 


D 


1  + 


JLLe. 


Y  ♦  { 


1  +  a 

Or 


)y  i 


(34) 


Thus,  the  detonation  velocity  is  seen  to  be  a  function  of  v0  ,  the 
initial  specific  volume,  corresponding  to  the  well-known  experimental 
result  for  solid  explosives. 
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EVALUATION  OP  PARAMETERS 

The  three  undetermined  parameters,  Y,  a,  and  b/Q,  which  appear  in 
the  LSZK  equation  of  state,  must  be  evaluated  by  using  experimental  data. 

It  can  be  seen  from  equation  (14),  which  describes  the  isentropic  P-p 
relation  for  an  LSZK  substance,  that  if  J-  *  or  <  y ,  p  ^  \  approaches 

1  ♦  Of  * 

Kp  a  as  p  approaches  zero.  We  assume  that  in  the  limit  of  low  pres¬ 
sures  the  detonation  products  behave  as  ideal  gases,  with  a  constant 
value  of  the  specific  heat  ratio,  denoted  here  by  *.  (A  reasonable  value 
for  H  seems  to  be  1.3^.  obtainable  by  averaging  the  ganmas  for  the  various 
gaseous  constituents  according  to  the  composition  of  the  products  at  low 
pressure.)  It  is  thus  clear  that  in  order  to  obtain  the  correct  behavior 
of  the  detonation  products  at  low  pressures,  we  must  set  a  **  1  . 

*-l 

The  remaining  constants  may  be  evaluated  by  referring  tc  the  experi¬ 
mental  results  for  the  dependence  of  the  detonation  velocity  on  the 
density.  After  a  particular  value  is  assigned  for  Y  (Y>* ;  ,  a  series  of 
values  for  b/Q  may  be  obtained  by  carrying  out  a  point  by  point  comparison 

1  B 

of  ^he  theoretical  plot  of  In  D  vs.  (lnq,  ♦  — In  — )  (obtainable  from 

equations  (31)  end  (3M)  with  the  experimental  plot  of  In  D  vs.  In  q, . 

Since  B/Q  must  be  a  constant,  the  accuracy  of  the  fit  is  determined  Dy 
the  amount  of  variation  in  the  values  of  B/Q  obtained,  and  Y  may  be  adjusted 
to  make  this  variation  a  minimum. 

This  process  has  been  carried  out  for  TNT,  using  an  empirical  relation 
between  detonation  velocity  and  loading  density  determined  at  the 
Explosives  Research  Laboratory,  at  Bruce ton. 8  This  relation  may  be  written 
D  ■  O.IT85  ♦  0.3225  >  where  D  is  in  centimeters  per  microsecond  and  ^ 

is  in  grams  per  cubic  centimeter.  Using  *  ■  1.3^*  and  with  the  heat  of 
detonation  chosen11  to  be  10l8  cal/gm,  the  results  are  V  »  2.78,  B/Q  ■ 
O.53562  and  o  -  2.9412. 

Equations  (31)  and  (34)  may  now  be  utilized  to  provide  the  dependence 
of  the  detonation  velocity  on  the  loading  density,  by  letting  the  parameter 
y  run  through  a  rarge  of  values.  Table  1  gives  a  comparison  of  this 
theoretical  curve  with  the  enpirical  relationship,  and  it  can  be  seen  that 
the  fie  is  quite  good. 

It  is  interesting  to  note  that  this  formalism  predicts  an  upper 
density  limit  to  the  deton&bility ,  at  tht  loading  density  ^  -  1.793  gia/cc 
and  detonation  velocity  D  ■  0.757  cm/ usee.  This  comes  about  because  of 
the  fact  that  at  this  point  the  value  of  the  parameter  y  is  zero  (y 
decreases  with  increasing  loading  density),  and  y  cannot  be  negative, 
because  of  its  physical  meaning  as  a  ratio  of  pressures  (see  equation  (7))- 
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In  fact,  equation  (7)  lilies  that  at  this  limiting  point,  the  only 
contribution  to  the  pressure  comes  from  the  elastic  part,  while  the  thermal 
pressure  vanishes.  Zeldorich  and  Kompaneets8  have  the  following  to  say 
about  the  physical  significance  of  this  phenomenon:  "It  is  possible  to 
have  charge  densities  for  which  the  thermal  energy  is  much  smaller  than 
the  elastic  part.  This  corresponds  to  y  being  nearly  zero.... It  is  not 
quite  clear  what  happens  when  the  charge  density  is  large.  It  can  be 
assumed  that  in  this  case  the  dissociation  reaction  does  not  go  to 
completion,  since  the  supply  of  chemical  energy  is  insufficient  for  over¬ 
coming  the  work,  required  by  the  elastic  repulsion  forces  between  the 
molecules.  It  appears  as  though  the  chemical  energy  does  not  suffice 
for  the  molecular  rearrangement  which  leads  to  an  explosion."  This 
prediction  is  especially  interesting  inasmuch  as  it  is  known  that  TNT 
exhibits  increasing  resistance  to  detonation  with  Increasing  loading 
density,  as,  in  fact,  do  most  solid  explosives. 

Now  that  values  for  the  constants  in  the  LSZK  equation  have  been 
arrived  at,  the  conditions  at  the  Chapman- Jouguet  state  may  be  computed 
by  using  the  expressions  developed  in  the  preceding  section.  This  has 

been  done  for  TNT  at  several  loading  densities  and  the  results  are 
presented  in  Table  2.  (The  temperatures  were  calculated  using  Cv  "  0.3 
cal/gm-iegree . )  Table  3  presents  experimental  values  for  the  C-J  state, 
determined  by  Dremin,  et  al;*“  the  correspondence  between  calculated 
values  and  experiment  appears  good,  except  for  the  temperatures,  which 
seem  tc  be  far  below  the  generally  accepted  values  for  detonation  tem¬ 
peratures  of  several  thousand  degrees. 

The  fact  that  the  ratio  of  elastic  pressure  to  total  pressure 
increases  as  the  loading  density  increases  may  be  verified  in  the  last 
colunn  of  Table  2.  For  instance,  for  ■  1.625,  this  ratio  is  0-97^, 
which  means  that,  for  the  isentrope  given  in  equation  (14),  97-4^  of  the 
pressure  comes  from  the  elastic  pressure  term.  For  ■  1.00  gra/cc,  this 
ratio  is  O.805. 

Conseq’iently ,  in  the  vicinity  of  the  Chapman -Jouguet  state,  the  LSZK 
isentrope  may  be  approximated  by  a  polytropic  relation,  with  exponent 
equal  to  2.78.  In  this  connection,  it  is  interesting  to  note  the  experi¬ 
mental  results  of  Deal,8  which  indicate  that  the  explosion  products 
isentrope  for  RDX-TNT  may  be  fitted  quite  closely  to  a  poly tropic  P-p 
relation,  with  V  ■  2.77 >  at  least  down  to  ^00  bars.  It  thus  seems  likely 
that  the  LSZK  equation  of  state  for  TNT  not  only  yields  the  proper  D  vs 
relationship,  but  also  provides  the  proper  isentrope,  both  in  the 
vicinity  of  the  Chapman-Jeaguet  state  and  in  the  larg^  expansion  limit 
of  low  pressure  and  density- 
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FLOW  FIELD  BEHIND  DETONATION  SHOCK 

The  isentropic  flow  behind  a  detonation  shock  in  a  spherical 
explosive  is  governed  by  the  differential  equati ons : 4 


du 

2  uc3 

06) 

* 

bJ* 

t 

2 

buA 

1 

-  c^} 

dc3 

-  2uca(l 

*  ?!  r 

(37) 

-  c3  k  * 

wnere 


u  “  particle  velocity 


radial  distance  of  detonation  shock  from  origin 


ff  n  _  f 

3  t*  t 


time 


o  *  sound  speed 
d3? 


f  "  {dpfe}  57  '  vhere  p  *  denslty  ‘  v  • 
Calculating  f  by  means  of  (14),  we  obtain: 

-  +  BY  (Y  -  l)p 


Y  -  1  -  - 
'  a 


K(i-tjSL)  +" 

a 


BYp 


or 


(38) 


(39) 


To  utilize  (39)  in  the  system  of  differential  equations  (36),  (37), 
it  is  necessary  to  express  f  as  a  function  of  c3 .  This  may  be  done 
(in  principle)  by  solving  (17)  for  p  in  terms  of  c3 ,  and  substituting 
the  result  into  (39)-  Unfortunately,  it  is  not  possible  to  invert 
(17)  analytically ,  in  closed  form,  so  that  an  alternative  approach 
must  be  used.  The  procedure  chosen  here  is  to  convert  equations 
(36)  and  (37)  into  a  set  in  which  the  dependent  variables  are  u 
and  p,  rather  than  u  and  c*  .  In  this  case,  we  can  use  f  in  the 
form  (39),  as  a  function  of  p.  To  effect  this  change  of  variable 
we  make  use  of  the  equation 


dc3  m  dc3  dp 
d£  “  dp  d§ 


(40) 
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Since  c3  *  dp/dp  in  the  isertropic  flow,  we  my  put  dca/dp  ■  d3P/dp? 
and  using  (17)  we  get: 


§£■  -  |  (^-2)  P  *  +  BV  (V  -  l)pV  “  2  (41) 

up  O  H 

\  -  J 

Consequently,  ■(£  (- - — )  p  01  +  BY  (V  -  l)p^  2]  ,  and  the 

Of  Of 

differential  equations  become : 


1 

where  c 3  ■  K  (1— ^)p  a  ♦  BYpV  1,  and  f(p)  is  given  by  (39)*  These 

equations  are  to  be  solved  subject  to  the  conditions  u  *  vl,  p  ■  n 
at  l  ■  D,  where  D 

D  ■  detonation  velocity 

Up  *  particle  velocity  at  the  detonation  shock 

Pp  *  density  at  the  detonation  shock. 

Up  and  pp  may  be  found  from  (33)  and  (25),  after  y  lias  been  found  from 
equation  (31)- 

These  calculations  have  been  carried  out  for  TK?  on  an  IBM-7090 
electronic  computer,  for  the  loading  densities  1.625,  l-59»  1-45,  1*30, 
1.14,  and  1.00  g m/cc.  The  results  are  given  in  Tables  4-9-  The  first 
column  's  a  dimensionless  distance,  the  radius  of  the  original  charge 
being  taken  as  the  unit.  Pressures  are  in  megabars,  velocities  in  centi¬ 
meters  per  microsecond,  energy  densities  in  megabar-cc  per  gram  and 
densities  in  grams  per  cubic  centimeter.  It  will  be  seen  that  the  para¬ 
meters  vary  in  the  well-known  way  first  demonstrated  by  Taylor,*  with  the 
region  of  constant  state  surrounding  the  origin. 
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CONCLUDING  REMARKS 

Up  until  now,  very  little  has  been  said  about  the  temperature.  To 
evaluate  this  quantity,  we  must  know  the  value  of  Cy,  which  is  not 
determined  by  the  other  ccnstants.  (Only  the  ratio  Cy/Cy^  is  determined 
by  equation  (^).)  If  Cv  is  taken  to  be  0.3  cal/gm,  ar,  approximate  average 
value  for  detonation  products,  the  C-J  temperature  (for  q,  ■  I.625  gm/cc) 
turns  out  to  be  582. 9° K,  which  seems  to  be  too  low.  This  is  connected 
with  the  phenomenon  of  the  decreasing  importance  of  the  thermal  pressure 
with  increasing  loading  density,  which  was  mentioned  above.  Though  this 
phenomenon  is  consistent  with  the  known  resistance  to  detonation  of  TNT  at 
high  densities,  and  with  the  experimental  results  of  Deal,8  it  is  not  yet 
certain  whether  it  is  a  real  effect  or  whether  it  is  a  result  of  the 
incompleteness  of  the  13ZK  theory*  In  any  case,  it  is  believed  that 
in  all  applications  where  the  temperature  is  not  needed,  and  only  an 
(E,  p,  v)  equation  of  state  is  required  (such  as  the  calculation  of  the 
non-reactive,  isentropic  expansion  of  detonation  products  by  means  of 
hydrodynamic  computer  codes),  the  LSZK  equation  of  state  (in  particular, 
equation  (3))  nay  be  used  with  confidence. 

It  is  probably  not  possible  to  decide  on  the  correctness  of  the 
U3ZK  equation  of  state  by  experimental  observations  of  the  detonation 
process  alone.  A  possible  approach  is  to  use  the  results  for  the  distri¬ 
bution  behind  the  detonation  as  initial  conditions  for  a  hydrodynamic 
code  computation  of  the  detonation  of  a  sphere  of  TNT  in  air,  using  the 
D5ZK  equation  of  state  for  the  expanding  detonation  products.  The  motion 
of  the  second  shock  through  the  product  gases  is  expected  to  be  a 
sensitive  function  of  the  equation  of  state  used,  and  one  can  attempt 
to  compare  the  calculated  results  with  the  evidence  obtained  from 
photographic  records.  The  behavior  of  the  air  shock,  though  a  much  less 
sensitive  function  of  the  equation  of  state  for  explosion  products, 
might  also  provide  a  useful  check. 

Preliminary  hydrodynamic  calculations  have  already  been  carried  out 
on  an  IBM-7090  and  are  reported  elsewhere;'0  more  refined  computations 
are  in  progress  at  the  present  time. 


*  Jacobs13  has  pointed  ou'  that  a  re mterpretat .on  of  the  partition 
between  elastic  and  thermal  energy  leads  to  a  theory  wh.ch  does  r,ot 
involve  a  limiting  density,  or  vanishing  thermal  pressure..-  This 
theory  retains  the  LSZK  form  for  the  equation  of  state,  but  does 
not  make  use  c *"  Zeldovich's®  arguments  for  the  physical  meaning  of 
the  cot  slants  ,  and  Cy. 
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Table  1 

Comparison  of  Detonation  Velocities  Calculated  for 
LSZK  Substance  with  Detonation  Velocities  Determined  at  Bruce ton 


*<£» 

D(—  );  ISZK 
'usee 

D(  -C-m  );  (Bruci 
'usee' 

1-7935 

0.7572 

0.7569 

1.6620 

0.7146 

0.7145 

1-5535 

0.6795 

0.6795 

1.4412 

0.6433 

0.6433 

1.3655 

0 . 6189 

0.5189 

1.2995 

0.5977 

0.5976 

1.2412 

0.5791 

O.5788 

1.1773 

0.5568 

0.5582 

1.1320 

0.5444 

O.5436 

1.1009 

0.5345 

0.5335 

1.0034 

0.5039 

0.5021 

0.9590 

0.4900 

0.4fi78 

O.9256 

0.4797 

0.4770 

0.9010 

0.4720 

0.4691 

O.8565 

0.4584 

0.4547 

0.8082 

0.4437 

0.4391 

0.7703 

0.4322 

0.4269 

0.7331 

0.4211 

0.4149 
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Table  2 


Detonation  Parameters  Calculated  with 
LSZK  Equation  of  State ,  for  TNT,  Q  ■  :>0l8  cal/gm 


{^(gm/cc) 

P(Kbars) 

[  mega bar -cc\( 
gram  ' 

p<g> 

u(-SH-) 

vusec 

d(-2L) 

'usee 

/degree^ 

1 'Kelvin ' 

1.625 

214.3 

0.06022 

0.188 

0.703 

582.9 

1-59 

203.5 

0.05973 

O.185 

0.691 

698.4 

I.45 

163.8 

0.0579 

1.986 

0.175 

0 . 646 

1141.7 

1-30 

127.6 

0.05607 

1.792 

0.164 

0.598 

1582.5 

1.14 

95-4 

0.05431 

1.583 

0.153 

0.547 

2013.3 

1.00 

72.2 

0.05293 

1.400 

0.144 

0.503 

2356.7 

1 
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Table  3 

Experinental  Values  for  Detonation  Parameters 
of  TUT  (Dremin,  et  al) 


*<«) 

'cc 

P(kbars) 

t  cm  x 

u( - ) 

'usee' 

D(— ) 
u*ec' 

1-59 

202 

O.I83 

0.180 

0.694 

1.45 

162 

0.172 

0.168 

0.650 

I.30 

123 

O.I56 

O.I56 

0.600 

1.14 

0  0 

Si 

0.557 

1.00 

64 

0.130 

0.129 

0.510 
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Table  4 

Detonation  Wave  for  TNT  (q>  ■  1. 

DISTANCE  VELOCITY  CENSITY 

X/RAOIIS  C8/USEC  GP/CC 


0. 

4.57265E-02 

9.  14531E-32 

1.37180E-C1 

1 . 82906E- 91 

2.28633E--0I 

2.74359E-01 

3.20086E-C1 

3.65812F-01 

4.11539E-01 

4.57265E-01 

4.8573CE-01 

5.14195E-01 

5. 42659E-01 

5. 71124E-01 

5.99589E-01 

6.28053E-91 

6.56518E-01 

6*8498  3E-01 

7.1 3447E-01 

7*4191 2E-01 

7 . 70377E-91 

7 .9884  1E-C1 

8.27396E-01 

8.55771E-H 

8.84235E-' 1 

9.  C0523E-01 

9.154776-01 

9, 291 05E-01 

9.4  14  17E-01 

9.52435E-01 

9.62185E-01 

9. 7070  IE-01 

9.78C22E-01 

9.84193E-01 

9. 89262E-01 

9.93282E-01 

9. 96309E-91 

9.98398E-01 

9.99609E-01 

l.COOOOE  00 


C. 

C. 

C. 

r 

V*  « 

c. 

c. 

c. 

c. 

c. 

c. 

5. 62295E— 05 
2.864C9E-&3 
6. 72  565E-G3 
1*11 8116-02 
1.60719E-02 
2.13 167L-02 
2.68  7246 -02 
3.27201E-02 
3. 88586E-02 
4.53  024E-0? 
5.20822C-02 
5. 92466E-02 
6.68685E-02 
7.50551E-02 
8. 39 o 76 E-02 
9.  38oC4C-02 
1 • 0)1 1  BE -01 
1. 06 3 75 C -01 
1. 126326-01 
1.1R890C-01 
1.25 147E-01 
1 . 31 4C5E-01 
1.37  662E -0 1 
1.43  919E-01 
1. 501774-01 
1.56434E-01 
1.62691E-01 
1.68 949E-01 
1. 75206C-01 
1.81463E-01 
1  •  87  72 1 E  -0 1 


1.29163E  Of 
1  • 29 lo 3 E  00 
i •  29 1 63E  OC 
1.29163E  0" 
1.29163E  0: 
1.29163E  0,* 
1 .29163c  00 
1.29163E  00 
1.29163E  00 
1.29163E  OC 
1.29163E  00 
I  •  30 326 E  OC 
1.3 1989 E  OC 
1.33975E  Ou 
1. 36220 E  OC 
1.3 8690 E  0. 
1.41 3t>7t‘  00 
1.44242E  Of 
1 .4731 1 E  00 
1.50582E  00 
1.540O7E  00 
1 . 5 7788 E  Oj 
1.6177 7 E  OC 
1. 66085 E  OC 
1.707b5E  0^ 

1 .75995  C  Of 
1 .792  78  £  Of 
1.82544E  Of 
1.85790 C  00 
1 .890  10 E  o: 
1.9  2200 E  00 
1.95337E  00 
1.98476  E  Of 
2 .0 1555 E  OC 
2.04590E  00 
2 .07578 E  Of 
2.10517E  Of 
2.1 3405  E  00 
2. 16 240 E  00 
2 . 19020 E  00 
2 . 2 1 743  E  Of 


d?5  go/cc) 


PHESSU^l  energy  oens 

MEGABAKS  f  r EG-CC ) /GM 


4.91537E- 02 
4.915  37E-02 
4.91537E-02 
4.91537E-C2 
4 .9 1 537E-C2 
4.9  1537E-02 
4.91537F-02 
4.91537E-C2 
4.91537E-02 
4 .9 1537E-02 
4.91537E-02 
5.C35R8E-02 
5.21 15  3E-02 
5.42632E-02 
5 . 6  7594E-C2 
5.959C 5 E-02 
6.2759 16-02 
6.62796E-02 
7 .C 1775E-02 
7.44899C-02 
7.92682E-02 
6 .458286-02 
9 .05  30  3  E - 0  2 
9.72468E-C2 
1.04932G-01 
1  .  13896E-01 
1  .  19769E-01 
1 .25843E-01 
1.32G5.E-Q1 
1 .38399E-C1 
1 .44880E-01 

1.5  148  2E-01 
1.58  192E-01 
1 .6500 JE-01 
1.71892E-01 
1.78856E-01 
1 .8588 3C-01 
1 .929526-01 
2.00060E-01 
2 .07 19  IE-01 
2 . 14333C-G1 


2.63C62E— C2 

2.63062F-C2 

2.63062E-C2 

2.63r62F-C2 

2. 63 C 62 F- 02 

2.63C62C-C2 

2.63C62E-32 

2.630  62p- 02 

2.63' 62E-C2 

2.630626-02 

2.6  30  62E-C2 

2. 66499E-C2 

2.71453E-C2 

2. 77425r-C2 

2.84252F-C2 

2.91858E-C2 

3.0C2C8E-C2 

3.C9299F-C2 

3. 1S152E-C2 

3.25814E-C2 

3.41357F-C2 

3.53P90E-C2 

3.67566E-C2 

3.8261CE-C2 

3. 99349F-C2 

4. 183C1F-C2 

4. 30454E-C2 

4.42709E-C2 

4. 55C45F-02 

4.6  7442F-C2 

4.758  776- C2 

4.92332E-02 

5. C4  786E-02 

5.  1  722Ct.  -C2 
5. 256 15E-C2 
5. 41953E-C2 
5.54217F-02 
5. 66392E-C2 
5. 78461E-C2 
5. 9C4 1 IE-02 

6. C2228F-C2 
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Detonation  Wave  for  TNT  (og,  ■  1.59  gm/cc) 


DISTANCE 

X/RAOILS 


0. 

4. 560640-02 

9. 12 129E-C2 

1 • 368 1 9E - ?1 

1.82426E-01 

2 • 28032E-01 

2.73639E-01 

3. 19245E-C1 

3. 64851E-01 

4.10458E-01 

4. 56064E-C1 

4.84596E-01 

5.  13127E-CI 

5.41659E-01 

5.70I9IE-CI 

5.98722E-01 

6.27254E-C1 

6.55785E-C1 

6.84317E-"1 

7.  12843E-01 

7.4 1 38CE-01 

7.69912E-G1 

7.98443E-C1 

8.26975E-01 

8. 55506E-C1 

8. 84Q38E-01 

9  .  C  0  3  4  8  C  -  0 1 

9.I5325E-01 

9.28974E-01 

9.4  1307E-01 

9.52344E-01 

9.62112E-0I 

9. 70643E-0 1 

9.77978E-01 

9.84161E-01 

9.8924CE-01 

9. 93269E -01 

9. 9630 IE-91 

9.98395E-0I 

9.99609E-CI 

I.C0000E  0C 


VELOCITY 

CR/USEC 


C. 

r» 

L  • 

c. 


c  • 

C  • 

c. 
c  • 
c. 

9.  91  348t  -05 
2.72H87E-03 
6.62  p94t-03 
1.09217E-02 
1.57492E-02 
2.09282E-02 
2.641546-02 
3.21  91 4  E -02 
3.82  54 6 E -02 
4.46 192E-Q? 
5.13 151C-02 
5.R3903L-C2 
6.59166E-02 
7. 39999E-02 

H. 27990t-02 
9.  25 653E-02 
9.  87  363E-C2 

I. 04  9  C7C-01 
1.11C78E-01 
1.17249E-01 
1.2342OE-01 
1.29591E-01 
1.35  762C-01 
1.41 933L-C1 
1.48 104E-01 
1.54275E-01 
1.60446L-01 
1.66617E-01 
1.72789E-01 
1.7R960C-01 
1.85131C-01 


CENS  I T Y 

GH/CC 

1 .262  7  1 E 

C. 

1.26271E 

0. 

1.26271E 

ov 

1 .26271 E 

Gu 

1.26271E 

Gu 

i .2627  1 E 

00 

1.26271E 

OC 

1 .2627  IE 

Gl 

1.26271E 

00 

1.26271 E 

00 

1 . 26271 E 

00 

1.27354E 

Ou 

1.28983E 

o: 

1.30925E 

0«J 

1.33128E 

00 

1.35553E 

00 

1.38182E 

oc 

1.41005E 

00 

1 .4402 1 E 

oc 

1.47235E 

Oc 

1.50660E 

00 

1.54317E 

oc 

1.58238E 

00 

1.62471E 

00 

1  .6709  l  E 

00 

1*7221 3E 

00 

1.7543  IE 

oc 

1.78645E 

Oc 

1.8  1 8  3.7  E 

00 

1.84995E 

oc 

1.881296 

Ou 

1.9123)E 

00 

1 .94294  E 

00 

1 .9 7  319 E 

00 

2  .90300  E 

00 

2 .0  3236  E 

00 

2  .06  124  E 

oc 

2 .0896 1 E 

oc 

2.11746E 

Ou 

2.14477E 

00 

2.17153E 

OC 

PRESSURE 

MECABARS 


4.67753E-C2 
4.67753E-02 

4 .6775  3E-02 
4.67753E-02 
4.67753E-02 
4  .677536-02 

4 .6775  3E-C2 
4.67753E-C2 
4.67753E-02 

4 .6775  36-02 
4.67753E-02 
4.78596E-02 
4.95 179E-02 
5 . 15509E-02 
5.39  16  2E-02 
5.66G35E-0 2 
5.96061E-02 
6.29463E-02 
6.66453E-02 
7.07384E-02 
7.52742E-02 
8.C3196E-02 
8 . 5966 3E-02 
9.2343RE-02 
9.96419E-02 
1 .08  155E-01 
1  •  13746E-G1 
1.19491E-01 
1.25382E-01 
1.31498 E- 01 
1 .375596-01 
1 .43825E-01 
1.50 195E-01 
1.56656E-01 
1.63 199E-01 
1 .698  10E-01 
1.76479E-01 
1.83193E-C1 
1 .8994  IE-01 
1 .96712E-01 

2.G3493E-01 


ENERGY  PENS 
(REG-CO/GM 


2 . 671 C8E-C2 
2.671C8F-C 2 
2.671C8E-C2 
2.671C8E-02 
2.671C8E-G 2 
2.671C8E-C2 
2.671CRE-C2 
2.671C8E-C2 
2.671CPE-C2 
2.671C8E-C2 
2.671 C8C-02 
2.7C2S4E-C2 
2. 7 5111 E- 02 
2. 8C933E-02 
2.87595E-C2 
2.95C19E-02 

3.C3171E-C2 
3.  12047F-0 2 
3. 2 1667E-C2 
3.32074E-02 
3.4334CE-C2 
3. 55670E-02 
3. 68914E-02 
3. 835876-C2 
3.99912E-C 2 
4. 1839CE-C2 
4. 3C226E-C2 
4.4216CE-02 
4. 541 71E-02 
4. 6624 IF-02 
4. 78347E-C2 
4.90471E-C2 
5.02593E-C 2 
5.  14694E-C2 
5. 26756E-C2 
5. 38762E-C2 
5.5C696F-02 
5. 6254CE-02 
5. 74282E-02 
5. 859C6E- 02 
5.97400E-G2 


warn  64-ao 


wo*  6 

Drton*tlo©  Wbv»  tor  TKt  (fe  *  1.45  w/cc) 


01  STANCE 
X/RA0IUS 


0. 

4.53543E-02 
9.  G7087E-02 
1.360636-01 
I.81417E-01 
2.267726-01 
2.721266-01 
3.17480E-01 
3.^28356-01 
4.081896-01 
4.535436-01 
4*821846-01 
5*108256-01 
5.394656-01 
5*681066^01 
5.96747E-01 
6*25387E-01 
5.54028E-01 
6*  82669E-01 
7*113096-01 
7*399506-01 
7*685906-01 
7.972316-01 
8.258726-01 
8.545126-01 
8*831536-01 
8*99565E-01 
9.146416-01 
9*283866-01 
9,408116-01 
9.519346-01 
9.617806-01 
9.703836-01 
9.777816-01 
9. 84018E-01 
9.891426-01 
9. 932076-01 
9.962676-01 
9.983806-01 
9.996056-01 
1.00000E  00 


VELOCITY 
CM/US EC 


0. 

0. 

c. 

0. 

0* 

0. 

0. 

0. 

0. 

1.71616C-04 
2.469*46-03 
6.055886-03 
1*02 157C-02 
1. 479066-02 
1.970006-02 
2.4901  IE -02 
3.03746L-02 
3.61180E-02 
4.21442E-02 
4.848G8 t-02 
5.517256-02 
6.228646-02 
6.99216E-02 
7.82267E-0 2 
8.743666-02 
9.32657E-02 
9.90948E-02 
1*049246-01 
1.107536-01 
1.166826-01 
1.22411E-01 
1.282406-01 
1.340696-01 
1.39899E-01 
1.45728E-01 
1.51 5576-01 
1.573866-01 
1.632156-01 
1,690446-01 
1.748736-01 


CENSITY 

GH/CC 

1.146616 

O'. 

1.146ol6 

0C 

1.146o IE 

00 

1.146616 

91 

l ■ 146olE 

00 

1 .146616 

o: 

1.146616 

1.146616 

*u- 

1.14661E 

00 

1.14661E 

00 

1.146616 

03 

1.155836 

cc 

1.17077E 

00 

1.188716 

Ou 

1.209G6E 

o: 

1  •  2  3 1 49  F 

00 

1 »255bl6 ■ 

OC 

1.28194E 

00 

1.30986F 

OC 

1.339&2E 

* 

1.37133E 

OC 

1.405196 

00 

1.441496 

p,  ** 

1.48069C 

OC 

1.52345E 

Ow 

1.570836 

O'. 

1.60C73F 

00 

1.630486 

OC 

1.66006F 

00 

1.689416 

OC 

1.71850F 

00 

1.747286 

OC 

1.77574E 

OC 

1.803836 

OC 

1.831526 

00 

1.85879E 

00 

1.885616 

OC 

1.911976 

oc 

1.937846 

00 

1.96321F 

rr. 

1.988C6E 

Oj 

PRESSURE 

KlGABARS 


3.BG136E-C? 
3.801366- 02 
3.801 36 E~ 02 
3.80136F-C2 
3.80136E-C2 
3.80 136E- 0.2 
3.801366-02 
3.8C 1 36 E- 02 

3. RG136E-D2 
3.80136C-C2 
3.00136E-C2 
3.68149E-C2 

4. vl36f!E-C2 
4.17629E-02 
4. 36 57 tC- 02 
4. 58094 E- 02 
4.821966-02 
5 .089896-02 
5.38663E-0? 
5.714986-0? 
6.078^36-02 
fe.48353F-02 
6. 9  364  "6-1.2 
7.44777 6-02 
8.03280F-C2 
8.714976-0? 
9. 163906-02 
9.625366-92 
1.009866- 01 
i.C  583'E-C'l 
1. i077bC-0i 
1.158146-01 
1.2G937C-CI 
1.261366-31 
1.314CCC-01 
1 • 3672  IE- 01 
1.420896-01 
1.47494E-01 
1.529266-01 
i.583‘r7E-01 
1.638366-01 


ENERGY  DENS 
(fcEG-CO/GH 


2  •  82422E-  r-2 
2.824226-02 
2 . 82*226-C  ? 
2 • 824226-C? 
2.82422E-02 
2.92422E-C2 
2.824  22E-0? 
2.824  226- r  2 
2 .824?2P-C  2 
2.82422E-0 2 
2.82422E-C2 
2 .85 }94E- ?2 
2.894526-02 
2 . 94732E-  >2 

3. CC777F-C2 
3.075136-02 
3. 14905E-02 
3.22947F- 0? 
3.31653F-C2 
3.41C63E-C? 
3.5 1237P-C? 
3.62268F-C? 
3.74287F-C2 
3. 874066-02 

4. C2149r-C2 
A.miPF-'*? 
4.?934?p-r2 
4.40050*'-  Z? 
4.5CB??F-C? 
4 .61642F- r2 
4.72491F-02 
4.83350F-T2 
4.94204F-:? 
6.C5" 34F-C? 
5. 15825F-C2 
5. 26662f-C? 
5.37229F-C2 
5.47P13F-C2 
5.5P3C  C6-C2 
5. 686796-02 
5.78936E-0? 


NOLTR  64-40 


Table  7 


Detonation  Wave  for  TNT  (po  ■  1-30  gVcc) 


DI STANCE 
X/RADI US 


0. 

4. 53858E-02 
9. 07715E-02 
L . 361 5  7E-01 
1. 81543E-01 
2.26929E-01 
2. 72315E-01 
3. 17700E-01 
3.63086E-01 
4. 084  72E-01 

4.  5385  8E-01 
4. 82401E-01 

5.  10944E-01 
5. 39487E-01 
5.68030E-01 
5. 96573E-01 
6. 251 16E-01 
6. 5365  9E-01 

6.  822  02E-01 
7. 10745E-01 
7.39288E-01 
7. 67831E-01 
7. 96374E-01 
8. 2491 7E-01 
8. 53460E-01 
8. 82003E-01 
8.98544E-01 
9.  13  74  8E-01 
9.27617E-01 
9.40161E-01 
9. 51395E-01 
9.61344E-01 
9.  7004  OE -01 
9.  7752  OE-Ol 
9.83828E-01 
9.89012E-01 
9. 93125E-01 
9.96222E-01 
9. 98361E-01 
9. 99600E-01 
1.00000E  00 


VELOCITY 

CM/USEC 


C. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

G. 

4.92819C-05 
2.52551E-03 
5.93218E-03 
9.86073E-03 
1.41695C-02 
1. 878 54 E -02 
2.36693E-02 
2.88034E-02 
3.41857E-02 
3.98  2  76 E -02 
4.57546E-02 
5.20076E-02 
5.86481E-02 
6.57666E-02 

7.  ''4  98  5 E  -0  2 
8.205691-02 
8.75274E-0 2 
9.29978E-02 
9.8468 3E -02 
1.03939E-01 
1  • 09 409 E -0 1 
1.  14880E-01 
1.20350E-01 
1.25821E-0 1 
1.31291E-01 
1.36762E-01 
1.42232E-01 
1.47702E-01 
1.53173E-01 
1.58643E-01 
1.64  114E-01 


CENS  ITY 
GM/CC 


1.0  21  38  E  00 
1.02138E  OC 
1.02138E  00 
1.02138E  00 
1 .02 1 38  E  00 
1.02138E  00 
1.02138E  00 
1 .02  1  38  E  00 
1 .0  2  l  38  E  00 
1.02138E  00 
1.02138F  00 
l .0  3098  E  OC 
1 .044  74  E  00 
1.06118E  00 
1.07977E  00 
1 . 100  22  E  00 
1.12239E  OC 
1.14619E  00 
1.17161E  00 
1.19870E  00 
1.22754E  00 
1.25833E  00 
1.291 33 E  00 
1.32694E  00 
1.36576E  00 
1.40873E  00 
t.  .4  36  1  3E  00 
1.46341E  00 
l  .49053  E  00 
1.51746E  00 
1.54416E  00 
1.57059E  00 
1 . 59672 E  00 
1.62251E  OC 
1.64795E  00 
1.67301E  00 
1.69765E  00 
1.72187E  OC 
1 .74565  E  00 
1.76896E  00 
l .79 180  E  00 


PRESSURL 

MEGABARS 


2.990C5E-02 

2 . 9900  5E-02 
2 .99005E-02 
2.99OO5E-02 

2 . 9900  5E-02 
2.99005E-02 
2.99005E-02 
2.99CG5E-02 

2 . 9900  5F-  02 

2 . 9900  5E-02 
2.99005E-C2 

3.C61C6E-02 
3.16466E-02 
3.29  137t-02 
3.43860E-02 
3.6G553E-02 
3. 79229E-C2 
3.99970E-02 
4 . 2292  IE-02 
4.48297E-02 
4.76  396E-02 
5.07624E-02 
5.42537E-02 
5.81922E-0? 
6.26922E-02 
6.7931 OE-02 
7.  1414CE-02 
7.4995PE-02 
7.86708E-02 
0.243  30E-02 
8.62759E-02 
9.01929E-02 
9.4l76oE-02 
9 t 82206E-02 
1.02317C-01 
1 .06457E-01 
1.  10636E-01 
1  .  1484  2E-01 
l.  19072E-01 
1.23316E-01 
1 .27567E-01 


ENERGY  DENS 
(MEG-CC ) /GM 


2.97234E-C2 
2. 97234E-C2 
2.97234F-C2 
2.97234F-C2 
2.97234E-C2 
2.97234E-C? 
2.97234E-C2 
2. 97234E-C 2 
2.97234E-C2 
2.97234E-02 
2.97234E-C2 
2. 99994E-C? 
3.03970E-C2 
3.08755F-C2 
3. 14212C-C2 
3. 2C2  76r- 02 
3.26914E-C2 
3. 34120F-02 
3.419C5F-C2 
3.5C3C2E-C2 
3. 59362E-C2 
3.69165E-C2 
3. 79822F-C2 
3.9  1499F- 02 

4.C4437E-C2 

4.  19C14E-02 
4.28445C-C2 
4.37944F-02 
4 . 4  74  96T- 02 
4. 57084E-C2 
4. 66693F- 02 
4.7  6306F- 02 
4. 85909F-C2 
4.95486E-02 
5.05024C-02 

5.  145C8F-C2 
5.23926E-02 
5. 33266E-02 
5.42515E-C2 
5. 51663E-C2 
5.6C70CE-C2 


o  >»■ 


NOLTR  64-40 


liable  8 


Detonation  Wave  for  TNT  (ft»  *1-14  gmjcc) 


Dl STANCE 
X/RADIUS 


• 

. 52631E-02 
9.C5262E-02 
1. 35789E-01 
1. 81052E-01 
2.26316E-01 
2. 71579E-01 
3. 16P42E-01 
3-621 05E-01 
4.07368E-01 
4. 52631E-01 
4.95417E-01 
5. 38202E-01 
5.80988E-01 
6.23774E-01 
6. 66560E-01 
7 . C9345E-01 
7.52131E-01 
7.94917E-01 
8. 37703E-01 
8.80488E-01 
8.93143E-01 
9. 05Q47E-01 
9.16197E-01 
9.26595E-01 
9. 36243E-01 
9. 45147E-01 
9. 5331 6E-01 
9.60760E-01 
9.67491E-01 
9#  73526E-01 
9. 78880E-01 
9.83573E-01 
9. 87625E-01 
9.91058E-01 
9.93896E-01 
9.96161E-01 
9. 97879E-01 
9.99074E-01 
9* 99773E-01 
i.OOOOOE  00 


VELOCITY 

CM/USEC 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

8.34041E-05 
3.96  314C-03 
9.32355E-03 
1.54886E-02 
2.22482E-02 
2.9526EE-02 
3.73264E-02 
4.571 34C-Q2 
5.48  296C-02 
6.49320E-02 
7.64985E-02 
8.03235E-02 
8.41 484E-02 
8.79733C-02 
9.17982E-02 
9.56232E-02 
9. 94  48  1 E  -0  2 
1.03273E-01 
1.07  098  E -0 1 
1. 10923E-01 
1.14  748E-0 1 
1.18573E-01 
1.22398E-01 
1.26223E-01 
1 . 30047E-0 1 
1.33872C-01 
1.37697E-01 
1.41 522E-01 
1.45  347E-01 
1.49172E-01 
1.52997E-01 


LENS  ITY 
GM/CC 


8.88453E-01 
8 .88453E-01 
8.88453E-01 
8.88453E-01 
8.88453E-01 
8.88453E-01 
8.88453E-01 
8.88453E-01 
8.88453E-01 
8  •  8845  3  E-0 1 
8.88453E-01 
9.03023E-01 
9.243 14  E-0 1 
9.49978E-01 
9.79245E-01 
1.01182E  OC 
1.047&9E  Ou 
1.08713E  00 
1.13074E  00 
1.17961E  00 
1.23582E  00 
1.25438E  00 
1.27290E  OG 
1.29136E  00 
1.30975E  00 
1 . 32804 E  OC 
1.34622E  00 
1.36429E  OC 
1.38222E  00 
1 . 40000 E  OC 
1.41761E  Oo 
1.43506E  OC 
1.45233E  00 
1 .46940 E  00 
1.48627E  00 
1.50292E  00 
1.51936E  00 
1.53557E  00 
1.55155E  00 
1 . 56729  E  00 
1.58278E  OC 


PRESSURE 

MEGABARS 


2.26994E-02 
2.26994E-02 
2.26994E-C2 
2 .26994E-02 
2.26994E-02 
2.26994E-C2 
2 • 26994T-02 
2.26994E-02 
2.26994E-02 
2.26994E-C2 
2.26994E-02 
2 . 35974E- 02 
2.49510E-02 
2.66488E-C2 
2 .8676QE-02 
3. 1C453E-02 
3.38000F-02 
3 . 7007  7E- 02 
4.07783E-02 
4 . 52926E-02 
5.C873  It-  02 
5.28096E-02 
5. 4789CE-02 
5.68U95E-02 
5.P8694E-02 
6 . 09666F- 02 
6 . 3099 2E-02 
6.5265  IE-02 
6. 7462  1C-02 
6.968/9E-02 
7.  19402E-02 
7.42167E-02 
7.65 149E- 02 
7.88  324E-02 
8.1 1668E-02 
8.35 1 56E-02 
8 • 5876  4E-02 
8.82467E-02 
9.0624CE-02 
9  •  3006  IE-02 
9  •  5390  5E-02 


ENERGY  DENS 
(MEG-CO/GM 


3.  11573E-02 
3.  11573E-02 
3. 11573E-C2 
3.1 1 573F-C2 
3.1 1573F-02 
3.1 1573E-C2 
3.1 1573E-02 
3.11 573E-C2 
3.1 1573F-C2 

3.  1 1573E- 02 
3. 11573E-C2 
3. 15776E-C2 
3.21966E-02 
3 . 29503E-C2 
3 . 382C0E-C2 
3.4P008E-C2 
3 . 5  P97CE- 02 
3.71218E-C2 
3. 84999E-02 
4 . 0C747F- C2 
4. 19256C-02 
4.25463E-02 
4. 317C2E-02 
4. 37967E-02 
4.44253E-C2 
4. 5C553E-C2 
4 • 56862E- C2 
4.631 73E-C2 
4 . 694R2F- 02 

4.  75782F- C2 
4. 82068E-02 
4.8P334E-02 
4.945  77E- 02 
5 . 00789E- 02 
5. 0696PF-02 
5. 13107E-02 
5.19203F-02 
5.25252F-G2 
5.3  1249E- 02 
5. 37191E-C2 
5. 43C73E- C2 


N01/FR  64-40 


Table  9 

Detonation  Wave  for  T!VT  •  1.00  gm/cc ) 


DI STANCE 
X /RADI  US 


0. 

4.46499E-02 

8.92817E-?? 

1. 33923E-01 

1 • 78563E-C 1 

2.23204E-01 

2.67945E-01 

3.12486E-31 

3.57127E-01 

4.01768E-01 

4.49272E-01 

4.77912E-:i 

5 . 0555  IE-9 1 

5. 351 90E-Q1 

5.63830E-G1 

5.92469E-01 

6. 21139E-J1 

6.49748E-01 

6. 78387  E-01 

7.07027E-01 

7.35666E-91 

7.64305E-:>1 

7.92945E-01 

8.21584E-CI 

8.50223E-01 

8.78R63E-C1 

8.95739E-01 

9.U281E-DI 

9.  25483E-01 

9.38347E-01 

9.49885E-01 

9.60118E-01 

9.69Q72E-01 

9. 76782E-91 

9.83290E-01 

9. 8B643E-91 

9.92892E-01 

9.95093E-01 

9.983O4E-01 

9.99586E-01 

l.OOOOOE  30 


VELOCITY 

Cw/JSEC 


C . 
0. 
0  • 

c  • 

.*> 

*.  • 


1.28S9K-J4 
2.15273c-03 
5. 162120-03 
8.  63  988  c-G  3 
1.24547t-0 2 
1.65385E-?2 
2.065420-02 
2.53841E-02 
3.01 246E-C  2 
3.53840C-02 
4.02828L-02 
4.57547l-C2 
5. 15 j05£-C2 
5.77454£-0? 
6.44  517E-02 
7.13  447  L-02 
7.66  343E-C2 
8.14240E-02 
8.6? 136E-C? 
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13  ABSTRACT 


Calculations  have  been  made  of  the  flow  field  in  the  isentropic 
region  behind  a  detonation  wave  in  TNT,  using  the  Landau-Stanyukovich 
equation  of  state  for  the  detonation  products  (as  described  by 
Zeldovlch  and  Kompaneets),  Adjustable  constants  in  this  equation 
have  been  evaluated  by  imposing  ideal  gas  behavior  on  the  detonation 
products  in  the  large  expansion  (low  density)  limit,  and  by  fitting 
to  an  experimental  curve  of  detonation  velocity  versus  loading 
density.  Calculated  values  of  Chapman-Joiguefc  variables  correspond 
fairly  well  with  experimental  values  at  various  loading  densities, 
with  the  exception  of  the  temperatures,  which  seem  to  be  far  too  low. 
This  is  connected  with  the  fact  that  the  theory  predicts  an  upper 
limit  to  the  loading  density  at  which  an  explosive  will  detonate;  at 
this  point  the  thermal  energy  vanishes  and  only  the  elastic  energy 
contributes  to  the  energy  of  detonation. 
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